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Abstract. An SIR epidemic model with free boundary is investigated. 
This model describes the transmission of diseases. The behavior of positive 
solutions to a reaction-diffusion system in a radially symmetric domain is 
investigated. The existence and uniqueness of the global solution are given 
by the contraction mapping theorem. Sufficient conditions for the disease 
vanishing or spreading are given. Our result shows that the disease will not 
spread to the whole area if the basic reproduction number Rq < 1 or the 
initial infected radius ho is sufficiently small even that Ro > 1. Moreover, 
we prove that the disease will spread to the whole area if Ro > 1 and the 
initial infected radius ho is suitably large. 
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1 Introduction 

Recently epidemic model has been received a great attention in mathematical 
ecology. To describe the development of an infectious disease, compartmental 
models have been given to separate a population into various classes based on 
the stages of infection |2|. The classical SIR model is described by partitioning 
the population into susceptible, infectious and recovered individuals, denoted by 
S, I and R, respectively. Assume that the disease incubation period is negligible 
so that each susceptible individual becomes infectious and later recovers with a 
permanently or temporarily acquired immunity, then the SIR model is governed 
by the following system of differential equations: 

S(t) = -pS{t)I(t)-inS(t) + b, 

I(t) = 0S(t)I(t)-fi 2 I(t)-aI(t), (1.1) 
R(t) = al(t) - fi 3 R(t), 

where the total population size has been normalized to one and the influx of the 
susceptible comes from a constant recruitment rate b. The death rate for the S, I 
and R class is, respectively, given by //i,/i2 an d ^3. Biologically, it is natural to 
assume that fi\ < min{/i 2 , ^3}. The standard incidence of disease is denoted by 
(3SI, where /3 is the constant effective contact rate, which is the average number 
of contacts of the infectious per unit time. The recovery rate of the infectious is 



denoted by a such that 1/oc is the mean time of infection. 



In 



13J, the threshold behavior was given. The authors showed that the basic 



reproduction number Rq 



bp 



determines whether the disease dies out 



n 



(Ro < 1) or remains endemic (Rq > 1). In [12] . a complete analysis of the global 
dynamics of an ordinary differential equation model with multiple infectious stages 
was presented, showing the same threshold behavior. For other works on various 
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13, 



20 



22 



26| 



types of SIR epidemic model, interested readers may refer to 
and the references therein. 

There are other compartmental combinations for modelling some other dis- 
eases. For example, the SI model describes a disease, such as herpes or HIV, with 
two stages, where individuals are infectious for life and never removed. The SIS 
model describes the case when individuals recover from the disease but there is no 
immunity, and they return to the susceptible class. Examples for this SIS model 
include sexually transmitted diseases, plague and meningitis. Unlike SIR models, 



SEI models 



191 ] assume that a susceptible individual first goes through a latent 
(exposed) period before becoming infectious. An example of this model is the 
transmission of SARS 2jJ, which is one of the serious diseases that human beings 
face at present. 

When the distribution of the distinct classes is in different spatial locations, the 
diffusion terms should be taken into consideration and thus an extended version 
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of the above SIR system (II. ip can be described as the following: 

' S t - d x AS = -f3S(t)I(t) - frStt) + b, x G n, t > 0, 

It - d 2 AI = pS(t)I(t) - fi 2 I(t) ~ al(t), x e n, t > 0, 

R t - d 3 AR = al(t) - fJ, 3 R(t), xett,t>0, (1.2) 

d v S = d v I = d v R = xe dVl, t > 0, 

S(x, 0) = S (x), I(x, 0) = I {x), R(x, 0) = R {x), x e H, 

where is a fixed and bounded domain in IR n with smooth boundary 90, and rj is 
the outward unit normal vector on the boundary. Here the homogeneous Neumann 
boundary condition implies that the above system is self-contained and there is 
no emigration across the boundary. The positive constants dj(z = 1,2,3) are the 
diffusion coefficients. 

It must be pointed out that the solution of system (ll.2p is always positive for 
any time t > no matter what the nonnegative nontrivial initial date is. It means 
that the disease spreads to the whole area immediately even when the infectious is 
confined to a small part of the area in the beginning. It doesn't match the observed 
fact that disease always spreads gradually. Recently the free boundary has been 
introduced in many areas, especially the well-known Stefan condition has been 
used to describe the spreading process. For example, it was used in describing the 
melting of ice in contact with water 



24j, in the mode 



6|, and in the dynamics of population 
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18 



21 



ing of oxygen in the muscle 



23| . There is a vast literature 



on the Stefan problem, and some recent and theoretically advanced results can be 
found in 4\. 

Motivated by the statements mentioned above, we are attempting to consider a 



SIR epidemic model with a free boundary, which describes the spreading frontier of 
the disease. For simplicity, we assume the environment is radially symmetric. We 
will investigate the behavior of the positive solution (S(r,t), I(r,t), R(r,t); h(t)) 

with r = |x| and i 6 K" in the following problem: 

f 3 t - ^AS = b- pS(r, t)I(r, t) - ^S(r, t), r > 0, t > 0, 

It - d 2 AI = pS(r, t)I(r, t) - // 2 /(r, t) - al(r, t), < r < h{t), t > 0, 

Rt - d 3 AR = al(r, t) - fi 3 R(r, t), < r < h(t), t > 0, 

< S r (0, t) = J r (0, t) = R r (0, t) = 0, t>0, (1.3) 

J(r, *) = R{r, t) = 0, r > h(t), t > 0, 

h'(t) = -fjil r (h{t),t), h(0) = h >0, t > 0, 

I S(r,0) = 5 (r), J(r,0) = J (r), R(r,0) = R (r), r > 0, 

where Aw = w rr + I] ^-w r , r = h(t) is the moving boundary to be determined, ho, di 
and \i are positive constants. The initial functions S , 1 and i? are nonnegative 
and satisfy 

( S e C 2 ([0, +oo)), J , R e C 2 ([0, ho]), 

\ J (r) = i2o(r) = 0, r e [/i , +oo) and J (r) > 0, r e [0, /i )- 

Ecologically, this model means that beyond the free boundary r = h(t), there is 

only susceptible, no infectious or recovered individuals. The equation governing 

the free boundary, h'(t) = —jil r {h{t),t), is a special case of the well-known Stefan 



condition, which has been established in 



2 If ] for the diffusive populations. 



The remainder of this paper is organized as follows. In the next section, we first 
apply a contraction mapping theorem to prove the global existence and uniqueness 
of the solution to the problem (II. 3p . Then we make use of the Hopf Lemma to 
give the monotonicity of the free boundary. Section 3 is devoted to prove that 



the disease will vanish if the basic reproduction number R < 1. In Section 4, we 
discuss the case Ro > 1. Our results show that for the case Ro > 1, the disease 
will spread to the whole area if ho is suitably large; while the disease will vanish if 
ho is sufficiently small. Our arguments are based on the comparison principle and 
the construction of appropriate supper solution of (11.31) . Finally, we give a brief 
discussion in Section 5. 

2 Existence and uniqueness 

In this section, we first prove the following local existence and uniqueness result 
by the contraction mapping theorem. We then use suitable estimates to show that 
the solution is defined for all t > 0. 

Theorem 2.1 For any given (S , 1 , R ) satisfying (II. 4p and any 7 G (0, 1), there 
is a T > such that problem ( II. 3p admits a unique bounded solution 

(S,I,R;h) G C 1+ ^ l+l)/2 {D^) x [C 1+7 ' (1+7)/2 (D T )] 2 x C 1+7/2 ([0, T]); 

moreover, 

||5'||c' 1 +^.( 1 +^)/ 2 (D|?) + H-^llc 1 +^'( 1 +^)/ 2 (D T ) + l|-^llc 1 +^'( 1 +^)/ 2 (D T ) + I l^llc 1+ ^/ 2 ([0,T]) — 

where = {(r,t) G M 2 : r G [0,+oo),t G [0,T]} and £> T = {(r,t) G M 2 : r G 
[0, /i(t)],tG [0,T]}. Here C and T only depend on h Q , 7, ||^o||c 2 ([o,oo)), ||^o||c 2 ([o,/i ]) 
and ||-Ro||c 2 ([o,/i ])- 
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Proof: We first straighten the free boundary as in 
C 3 [0,oo) satisfying 



3- 



Let £(s) be a function in 



£(s) = 1 if \s - ho\ < ^, £(s) = if \s - ho\ > ^, \e(s)\ < |- for all a. 
Consider the transformation 

(y, t) -)• (x, t), where x = y + £(|y|)(/i(t) " WM), V e 
which leads to the transformation 

(s,t) -)■ (r, t) with r = s + £(s)(h(t) - h ), < s < oo. 

As long as 

IM*)-J>ol< 

the above transformation x — >• y is a diffeomorphism from R n onto R n and the 
transformation s — > r is also a diffeomorphism from [0, +oo) onto [0, +oo). More- 
over, it changes the free boundary r = h(t) to the line s = h . Now, direct 
calculations show that 

ds 1 VMKt),s), 

B(h(t),s), 
C(h(t),s), 
D(h(t),s). 



dr l + ?(s)(h(t)-ho) 

8h Z"(s)(h(t)-h ) 

<9r 2 [i + £( s )(h(t)-ho)} 3 " 
1 ds 



h'(t)dt l + £'(s)(h(t)-h ) 
in - l)y/A 



s + Z(s)(h(t)-h ) 
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Now, if we set 



S(r,t) = S(s + Z(s)(h(t)-ho),t) := u(s,t), 

I(r, t) = I{s + £(s)(h(t) - ho),t) := v{s, t), 
R(r, t) = R(s + £(s){h{t) - h ),t) := w(s, t), 

then the free boundary problem (11. 3p becomes 

' u t - Ad 1 A s u - (Bdi + h'C + Ddi)u s = b- (5uv - /iiu, s > 0, t > 0, 

v t — /W 2 A s t> — (Bd2 + h'C + Dd2)v s = (3uv — n%v — av, < s < ho, t > 0, 

w t — Ad 3 A s w — (Bd 3 + h'C + Dd$)w s = av — Hzw, < s < ho, t > 0, 

u s (0,t) =v s (0,t) =w s (0,t) =0, t>0, 

v(s,t) = w(s,t) =0, s > ho, t > 0, 

h'(t) = -fiv s {h , t), h{0) = h , t> 0, 

k u(s, 0) = u (s), v(s, 0) = v (s), w(s, 0) = w (s), s < 0, 

where A = A(h(t),s), B = B(h(t),s), C = C(h(t),s), D = D(h(t),s) and u 



S ,v = Iq,w = Rq- 

We denote h* = -fiv' Q (h ), and for < T < 8(1 +V) > set 

H T = {h G C l [0,T] : h(0) = h , h'(0) = h*, \\ti - h*\\ C (io,T]) < l}, 

U T = |w G C([0,+oo) x [0,T]) : u(s, 0) = m (s), ||u - «o|U-»([o,+oo)x[o,T]) < 1 

Vt = [v e C([0,oo) x [0,T]) : = for s > h , < t < T, 

v(s,0) = w (s)for0 < s < h , \\v - Uo||l°°([o,oo)x[o,t]) ^ l}, 

W T = {w G C([0,oo) x [0,T]) : io(s,t) = for s > h , < t < T, 
w(s, 0) = w {s) for < s < h , \\w - wo||i/*>([o,oo)x[o,t]) < 1 



Noticing the fact that for hi, h 2 G H T , due to hi(0) = h 2 (0) = h , we have 



\\h - h 2 \\c([o,T}) < T\\h[ - h' 2 \\c([o,T]), (2.3) 

it is not difficult to see that IV := Ut x Vt x Wt x Ht is a complete metric space 
with the metric 

V((ui,vi, wi; hi), (u 2 ,v 2 ,w 2 ; h 2 )) = \\m - u 2 \\l^([o,+o Q )x[o,t}) 

+ \\Vi — f2||L°°([0,oo)x[0,T]) + \\Wl — W 2 \\l°°([0,oo)x[0,T]) + \\K — h 2 \\c([Q,T])- 

Next, we shall prove the existence and uniqueness result by using the contrac- 
tion mapping theorem. Applying standard L p theory and the Sobolev imbedding 



theorem 



151 ]. we can find that for any {u, v, w; h) G IV, the following initial bound- 



ary value problem 

Ut — Adi/S. s u — (Bdi + h'C + Dd\)u s = b — (3uv — niu, s > 0, t > 0, 

v t — Ad 2 A s v — (Bd 2 + h'C + Dd 2 )v s = f3uv — jj, 2 v — av, < s < h , t > 0, 

w t — Ad 3 A s w — (Bd 3 + h'C + Dd 3 )w s = av — ^w, < s < h , t > 0, 

u s (0,t) =v s (0,t) =w s (0,t) =0, t>0, (i ' U 

v(s,t) = w(s,t) =0, s > h , t > 0, 

k u(s, 0) = u {s), v(s, 0) = v {s), w(s, 0) = w (s), s < 
admits a unique solution 

(u,v,w) e ]C 1+ ^ 1+ ^ 2 ([0,+oo) x [0,T])] 3 



and 



|'"||c 1 +^.( 1 +'')/ 2 ([0,+oo)x[0,T]) ^ ^i, (2.5) 

P||c 1 +t.( 1 +t)/2([o,Ao]x[0,T]) ^ ^1) ( 2 -6) 

l*llc 1 +-r.( 1 +^)/ 2 ([o,Mx[o,T]) ^ ^1) (2-7) 
9 



where K x is a constant depending on 7, h , ||So||c*[o,+oo)j ll^ollc^o.M and 1 1 1 1 [o,Ao] - 
Now, we define h(t) by the sixth equation in (12.21) as the following: 



h(t) = ho - H I v s {h ,r)dT, (2. 



then we have h'{t) = —[iv s (h ,t), h(0) = h and h'(0) = —}iv' (ho) = h*. Hence 
h'(t) e C^ 2 {[0,T]) with 

\\h\t)\\ c ^ {m <K 2 ^ l iK l . (2.9) 

In what follows, we define a map 

T : I> — ► [C([0,+oo) x [0,T])] 3 x C\[0,T}) 

by J-"(m(s, t),v(s, t),w(s, i); h(t)) = (u(s, t),v(s, t), w(s, t); h(t)). It is obvious that 
(u(s,t),v(s,t),w(s,t):h(t)) G T T is a fixed point of T if and only if it solves (12.21) . 

Similarly as in 7|], there is a T > such that J 7 is a contraction map- 
ping in IV. It follows from the contraction mapping theorem that there is a 
(u(s, t), v(s, t), w(s, £); h(t)) in IV such that 

^(u^s, t),v(s, t),w(s, t); h(t)) = (u(s, t),v(s, t),w(s, t); h(t)). 

In other words, (u(s, t), v(s, t), w(s, t); h(t)) is the solution of the problem (12. 2p 
and thereby (S(r, t), J(r, t), R(r, t); h(t)) is the solution of the problem (II. 3p . More- 
over, by using the Schauder estimates, we have additional regularity of the solution, 
h(t) G C 1+ ^ 2 ([0,T]),S G C 2+ ^/ 2 ( (0, +00) x (0, T]) and J,i?G C 2+ ^ 1+ ^ 2 ((0,h(t))x 
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(0, T]). Thus (S(r,t), I(r,t), R(r,t); h(t)) is the classical solution of the problem 
tfH. □ 
To show that the local solution obtained in Theorem 2.1 can be extended to 
all t > 0, we need the following estimate. 

Lemma 2.2 Let (S, I, R; h) be a bounded solution to problem ( I1.3P defined for t G 
(0,T ) /or some T G (0, +oo]. T/ien there exist constants C\ and Ci independent 
ofT such that 

< S(r, t) < Ci for0<r < +oo, t G (0, T ). 

< J(r, t), i?(r, t) < C 2 for < r < t G (0, T ). 

Proof: It is easy to see that S > 0, / > and R > in [0, +oo) x [0, T ) as long 
as the solution exists. 

Using the strong maximum principle to the equations in [0, h(t)] x [0,To), we 
immediately obtain 

S(r, t), I{r, t), R{r, t) > for < r < h(t), < t < T . 

The upper bounds of the solution are followed from the maximum principle, 
we omit the proof here. □ 

The next lemma shows that the free boundary for problem (11.31) is strictly 
monotone increasing. 
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Lemma 2.3 Let (S,I,R;h) be a solution to problem (jl.3p defined for t G (0, To) 

/or some To G (0, +oo]. T/ien t/iere exists a constant C3 independent of To suc/i 

0<fc'(t)<C 3 /or fe(0,T o ). 

Proof: Using the Hopf Lemma to the equation of I yields that 

l r (h(t),t)<0 forO<t<T . 

Hence h'(t) > for t G (0,T ) from the Stefan condition. 

Next we show that h'(t) < C3 for all t G (0,T ) and some C3 independent of 
Tq. As in [2l| , we define 

tt = tt M := {{r, t) : h{t) - ikT 1 < r < /i(t), < t < T } 

and construct an auxiliary function 

w (r,t) := C 2 [2M{h{t) - r) - M 2 {h(t) - r) 2 ]. 

We will choose M so that w(r,t) > I(r,t) holds over Jl. 
Direct calculations show that, for (r, £) G f2, 

w t = 2C 2 Mh'(t)(l - M(h(t) - r)) > 0, 

-Aw = 2C 2 M 2 , PSI - (H2 + a) I < $C X C 2 , 

and then 

w t - d 2 Aw > 2d 2 C 2 M 2 > $C X C 2 in Q 
12 



if M 2 > On the other hand, we have 

w(h(t) - M~\t) = C 2 > I(h(t) - M~\t), w(h(t),t) = = I(h(t),t). 

Hence, if we can choose M such that I (r) < w(r, 0) for re [h — M _1 , ho], then we 
can apply the maximum principle to w — I over Q to deduce that I(r, t) < w(r, t) 
for (r, t) G Vl. It would then follow that 

I r (h(t),t) > w r (h{t),t) = -2MC 2 , h\t) = -fil r (h(t),t) < C 3 := 2MC 2 fi. 

To complete the proof, we only have to find some M independent of T such 
that I (r) < w(r, 0) for r £ [ho — M _1 , ho]. We calculate 

Wr (r, 0) = -2C 2 M[\ - M(h - r)] < -C 2 M for r e [ho - (2M)" 1 , h ]. 

Then upon choosing 




we have 



w r (r,0) < -MC 2 < 



h\W < J oW for re [ho- (2M)-\h ]. 



Since w(/io,0) = Io(h ) = 0, the above inequality implies that 



w(r, 0) > J (r) for re[/i - (2M)" 1 , /i ]. 



Moreover, for r G [/iq - M 1 , ho — (2M) 1 ], we have 



w(r,0) > -C 2 , J (r) < \\I \\c H[ oM) M < -f 2 . 



Therefore /o( r ) < ^( r ; 0) for r £ [h — M 1 , /i ]. This completes the proof. 



□ 
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Theorem 2.4 The solution of the problem ( II. 3p exists and is unique for all t G 
(0,oo). 

Proof: It follows from the uniqueness of the solution that there is a number T max 
such that [0,T max ) is the maximal time interval in which the solution exists. Now 
we prove that T max = oo by the contradiction argument. Assume that T max < oo. 
Then it follows from Lemma 2.2 that there exist C\,C<i and C 3 independent of 
Tmax such that for t G [0,T max ) and r G [0, h(t)], 

< S(r,t) < C 1} (r,t) G [0,+oo) x [0,T max ), 

< I(r,t), R{r,t) < C 2 , (r,t) G [0,h(t)} x [0,T max ), 

h < h{t) <h + C 3 t, < h\t) <C 3 ,te [0, T max ). 

We now fix Sq G (0,T max ) and M > T max . Then by the standard parabolic regu- 
larity, we can find C4 > depending only on 5 , M, C\, C2 and C3 such that 

C 1 +t[0,+oo) j 

for i G [<5o, T max ). It then follows from the proof of Theorem 2.1 that there exists a 
t > depending only on Cj(i = 1, 2, 3, 4) such that the solution of problem (II. 3p 
with initial time T maa: — r/2 can be extended uniquely to the time T ma;c — r/2 + r. 
But this contradicts the assumption and thereby the proof is complete. □ 

Remark 2.1 It follows from the uniqueness of the solution to (II. 3p and some stan- 
dard compactness arguments that the unique solution (S, I, R, h) depends continu- 
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ously on the parameters appearing in (jl.3p . This fact will be used in the following 
sections hereafter. 

We next decide when the transmission of diseases is spreading or vanishing. 
We need to divide our discussion into two cases: R < 1 and R > 1. 

3 The case R < 1 

It follows from Lemma 2.3 that r = h(t) is monotonic increasing and therefore 
there exists hoc G (0, +00] such that lim 4 _>. +00 h(t) = h^. The following theorem 
shows that the transmission of diseases is vanishing in the case that Rq < 1. 

Theorem 3.1 If R (=- < 1, then lim^ +00 ||/(-,*)||c([o,h(t)]) = and 

/ioo < 00. Moreover, lim t ^ +00 ||-R(-,£)||c([o,h(t)]) = and lim t _> +00 5(r,i) = ^ 
uniformly in any bounded subset of [0, 00). 

Proof: It follows from the comparison principle that S(r, t) < S(t) for r > and 
t G (0, +00), where 

S(t) := — + (N-S-olloa - — )e-^*, 
which is the solution of the problem 

^- = b- f x 1 S, t>0; 5(0) = 115-01 (3.1) 
at 

Since lim^oo S(t) = ^ , we deduce that 

limsup t _ >+tX) S(r,t) < ^ uniformly for r G [0, 00). 
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Recalling the condition R < 1, there exists T such S(r, t) < ^ 1 " L ^° in [0, oo) x 

[T , +00). Now I(r,t) satisfies 

I t -d 2 AI<[^^-^ 2 -a]I(r,t), 0<r<h{t), t > T , 
l(r,t) = 0, r = h{t),t>0, (3.2) 

I(r,T )>0, 0<r<h{T ). 

Therefore ||/(-, i)||c([o,/i(t)]) - > as t — > 00, since that j3^ — ^2 — « < 0. We 
then have \\R(-, t)\\c([o,h{t)]) — ► as t — > 00 from the third equation of fll.3p . 
Next we show that < +00. In fact, direct calculation yields 
d f m 



— / r n - l I(r,t)dr 



h{t) 

r n - l I t (r, t)dr + /i'(t)/i n - 1 (t)/(/i(t), t) 







r h(t) f h(t) 
/ d 2 r n_1 A/dr+ / I(r,t)(pS(r,t) - fi 2 - a)r n - l dr 
Jo Jo 

h(t) rh(t) 

d 2 (r n - l I r (r,t)) r dr+ / I(r, t)(/3S(r, t) - /i 2 - a^^dr 



jo 



A rh(t) 

-—h n - l h'{t)+ / I(r,t)([3S(r,t) - fi 2 - a)r n - l dr. 



V Jo 
Integrating from T to t (> T ) yields 

\ r n - l I(r,t)dr = / r™" 1 /^, T )dr + — h n (T ) - — h n (t) 
Jo Jo W n/j 

+ [ [ ° I(r,s)(pS(r,s)- fM 2 - ay^drds, t>T . (3.3) 

J T Jo 

Since < S(r,t) < for r e [0, h(t)) and t > T , we have 

PS(r, t) - fj, 2 - a < for t > T , 

hit) fh(T ) 1 > 

r"- 1 /(r, £)dr < / r"- 1 /(r, T )dr + — h n (T ) - —h n (t) for t > T , 



Jo nfi n/i 
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which in turn gives that < oo. 

Then it follows from the first equation of (11. 3p that lim^+oo S(r, t) = 
uniformly in any bounded subset of [0, oo). □ 

4 The case R > 1 

In order to study the case that the reproduction number R > 1, and for later 
applications, we need a comparison principle, which can be used to estimate S(r, t), 
I(r,t), R(r,t) and the free boundary r = h(t). As in |7j], the following comparison 
lemma can be obtained analogously. 

Lemma 4.1 Suppose that T G (0,oo), h G C\[0,T]), S G C([0,oo) x [0,T]) n 
C^i&oc) x (0,T\), 7, Re C(D* T )r\C 2 >\D* T ) with D* T = {(r,t) G M 2 : < r < 
h(t),0 <t<T}, and 

( S t - d x AS > b - i^tS, < r, < t < T, 

It - d 2 AI > {/3S -/i 2 - a)I, < r < h(t), < t < T, 

i? t - d 3 Ai? > aJ - a* 3 77, < r < fr(t), < t < T, 

5 r (0,t) > 0, 7 P (0,t) > 0, ^.(0,*) > 0, < t < T, 

7(r, t) = R(r, t) = 0, r > fc(t), < t < T, 

> -fil r {h{t),t), h(0) >h , <t <T, 

I S(r,0) > S {r), 7(r,0) > 7 (r), %, 0) > i? (r), < r < /i . 

TTien t/ie solution (S,I,R;h) of free boundary problem (II. 3p satisfies 
S(r, t) < S(r, t), h(t) < h(t) for r G (0, oo) and t G (0, T], 
7(r, t) < 7(r, t), R(r, t) < R{r, t) for r G (0, h(t)) and t G (0, T] . 
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Next we show that if h and \i are sufficiently small, the disease is vanishing 
for the case Rq> 1. 



d 
8M : 



, and 



Theorem 4.2 If R (:= > 1, h < min{^, ^} and p < 

then hoo < oo. Where ko = {3Ci — ^2 — a > 0, G\ = max 1 1 1 *So | |oo, ^-j 
M = |max{||/ ||oo, H-Rolloo}- 

Proof: We are going to construct a suitable upper solution to (II. 3p and then 



apply Lemma 4.1. As in 



7|, we define 



S(r,t) = C 1 , 



I 



R 



Me-^V{r/h(t)), < r < h(t), 

0, r>h(t), 

Me-^V(r/h(t)), < r < h(t), 

0, r>h(t), 



and 



h{t) = 2h {2 - e~*), t > 0; V{y) = 1 - y 2 , < y < 1, 
where Ci = max ^ 1 1 6*0 1 1 00 , j^j, 7 and M are positive constants to be chosen later. 
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Denoting k = fiC\— n 2 — a, we have k > since i? > 1. Direct computations 
yield 

S t - d x AS = > b - mS, 
I t -d 2 AI- {pS-H2-a)I 
= I t - d 2 AI - kj 

= Me-^l-^V -rhh V'-d 2 h V" - d 2 h V - k V] 

r 

R t - d 3 AR - (al - fi 3 R) > Me-*[ - -y - a] 

for all < r < h(t) and t > 0. On the other hand, we have //(t) = 2/i 7e~ 7 * 
and -nT r (h(t),t) = 2M/iF 1 (i)e^'. Moreover, it follows that 5(r, 0) > 5 (r), 
/(r,0) = M(l - ^) > |M, %,0) = M(l - ^) > |M for r G [0,/i ]- Noting 
that h(t) < Ah , we now choose 

4 

M = -max{||/ ||oo, H-Rolloo} 

and take 

d d 
7= 16M' ^~ 8M' 
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where d := min{di, d 2 } and h < min | yj j^, \J Then we have 

S t - d^S >b- pSL - fiiS, < r, t > 0, 

I t - d 2 AI > (pS -fjL 2 - a)T, < r < i > 0, 

Rt - d 3 AR >al- fi 3 R, < r < h(t), t > 0, 

S r (0,t) = 0, J r (0,t) > 0, i? r (0,t) > 0, t > 0, 

7(r, t) = i?(r, i) = 0, r > 7&(i), i > 0, 

h'(t) > -fil r (h(t),t), h(0) = 2h > h , t> 0, 
[ S(r,0) > S (r), /(r,0) > J (r), %, 0) > i? (r), r > 0. 

Hence we can apply Lemma 4.1 to conclude that h(t) < h(t) for t > 0. Therefore, 

we have /i^ < lim^oo = 4/i < oo. □ 

For the case that R > 1, we next prove that if /i is suitably large, the disease 

is spreading. 

Lemma 4.3 If < oo, then lim t _^ +00 ||/(-,£)||c([o,h(t)]) — 0. Moreover, we have 
lim t _^ +00 ||-R(-,t)||c*([o,A(t)]) = and lim t _^ +00 S(r,t) = — uniformly in any bounded 



subset of [0, oo). 

Proof: Assume limsup t ^ +oc ||J(-, i)||c([o,/i(t)]) = 5 > by contradiction. Then 
there exists a sequence (r k ,t k ) in [0, h(t)) x (0, oo) such that I(r k ,t k ) > 5/2 for all 
k G N, and — >■ oo as — >■ oo. Since that < < < /i^ < oo, we then have 
that a subsequence of {r n } converges to r G [0, /i<x>)- Without loss of generality, 
we assume r k — >• r as /c — )■ oo. 

Define S*(r,*) = S(r,£ fe + t), I k (r,t) = I(r,t k + t) and # fc (r,t) = R(r,t k + 
i) for r G (0,h(t k + t)),t G (— t k ,oo). It follows from the parabolic regularity 
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that {(S k ,I k ,R k )} has a subsequence {{S ki , I ki , R ki )} such that (S ki , I ki , R ki ) ->■ 

(S*, /, i?) as i — > oo and (5, 7, i?) satisfies 

S t - diAS = b- pSI - n ± S, < r < h^, te (-00, 00), 
I t - d 2 AI = /3SI - fx 2 I - eel, < r < h^, t E (-00, 00), 
R t — d 3 AR — al — 113R, < r < h^, t e (—00, 00). 

Since J(r , 0) > 5/2, we have J > in [0, h^) x (—00,00). Recalling that (f3S — 
(j,2 — ct) is bounded by M := /3max{^-, ||So||l°°} + A*2 + ol. Applying Hopf lemma 
to the equation I t — d 2 AI > —MI at the point (0, h^) yields that I^h^, 0) < — er 
for some cr > 0. 

On the other hand, is increasing and bounded. Moreover, for any < a < 
1, there exists a constant C, which depends on a, ho, \\Io\\c 1+a [o,ho] an< ^ hoo, sucn 
that 

||7|| C l+a,(l+a)/2([ 0i h(t)) x [0,oo)) + 1 1 ^ 1 1 C 1 +«/ 2 ([0,Oo)) < C- 

In fact, let us straighten the free boundary in a way different from that in Theorem 
2.1. Define 

then direct calculations yield that 

h '(t) r h o a t h 2 Q A 

therefore, v(s,t) satisfies 

vt ~ d 2 j^A s v - ^sv s = v(/3u - pL 2 - a), < s < h , t > 0, 
v s (0,t)=v(h ,t) = 0, t>0, (4.2) 

v(s,0) = I (s) > 0, 0<s</i . 
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m), we 



This transformation changes the free boundary r = h(t) to the fixed line s = h , 
at the expense of making the equation more complicated. It follows from Lemmas 
2.2 and 2.3 that 

h'(t) 

\\v(/3u - n 2 - a)\\ L ~ <M X , < M 3 . 

Applying standard L p theory and then the Sobolev imbedding theorem 
obtain that 

|p||c 1 +«.( 1 +«)/ 2 ([o,/i ]x[o,oo)) < C4, 

where C4 is a constant depending on a, h , Mi, M 2 , M 3 and ||/o||c 2 [o,/i ]- This ™" 
mediately leads to (14. ip . 

Since ||/i||ciW2(ro,oo)) — ^ an d ^(*) * s bounded, we then have h'(t) — > as t — >• 
00, that is, I r (h(tk),tk) — )• as — )• 00 by the free boundary condition. Moreover, 
it follows from the inequality ||/||c 1 + a .( 1 +«)/ 2 ([o,/i(t))x[o,oo)) < C that I r (h(t k ),t k + 0) = 
{Ik)r{h{tk),0) — )■ ^(ZiocO) as A; — >■ 00, which leads to a contradiction to the fact 
that 1^00,0) < -a < 0. Thus linu_> +00 ||i"(-,t)||c([Q ) A(t)]) = 0, and thereby 
lim t _> +00 ||i?(-,t)||c([o,fe(t)]) = and lim^+oo S(r,t) = ± uniformly in any bounded 
subset of [0, 00). □ 

Let \\(R) be the principal eigenvalue of the operator —A in Br (open ball with 
radius R) subject to homogeneous Dirichlet boundary condition. It is well-known 
that Ai(i?) is a strictly decreasing continuous function and 

lim \\{R) = +00 and lim XAR) = 0. 
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Theorem 4.4 If R (:= ,^ +Q s ) > 1, then = oo provided that h > h^, where 
Xi(K) = ^(R -l). 

Proof: Assume that hoo < +00 by contradiction. It follows from Lemma 4.3 

that lim t ^ +00 \\I{-,t)\\ C ([o,h(t)]) = °- Moreover, lim^ +00 S(r,t) = ^ uniformly 

in the bounded subset Bh . Therefore, for e > 0, there exists T* > such that 

S(r, t) > \ - e for t > T*, r G [0, h(t)). We then have that I(r, t) satisfies 

f J t - d 2 AI > -e)-^- a)I, < r < h , t > T*, 

I l r (0,t) = 0, I(ho,t)>0, t>T\ (4.3) 

[ I(r,T*) > 0, < r < h . 

It is easy to see that I(r,t) has a lower solution I_(r,t) satisfying 

It ~ d 2 A£ = 08(£ - e) - // 2 - a)/, < r < h , t > T\ 
J r (0,t) = 0, I(h ,t) = 0, t>T*, (4.4) 

l(r,T*) = I(r,T*), < r < h . 

Since /io > we can choose £ sufficiently small such that — e) — — a) > 

d2^i{ho), it follows from well-known result that / is unbounded in (0, ho) x [T*, 00), 

which leads to a contradiction that lim^+oo ||/(-,£)||c([o,/i(t)]) = 0. □ 

5 Discussion 

In this paper, we have considered the SIR epidemic model describing the transmis- 
sion of diseases and examined the dynamical behavior of the population (S, I, R) 
with spreading front r = h(t) determined by (jl.3p . We have obtained the asymp- 
totic behavior results. 
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The basic reproduction number R (= Ml (^ +Q! ) ) is important but not unique 
factor to determine whether the disease dies out or remains endemic. It is shown 
that if Rq < 1, vanishing always happens or the disease dies out (Theorem 3.1). 
If Rq > 1, spreading happens provided that h is sufficiently large (Theorem 4.4) 
and vanishing is possible provided that ho is small (Theorem 4.2). 

We feel it is reasonable to conclude that (jl.3p is promising alternatives to (jl.ip 



and (ll.2p for the modeling of disease spreading, and there is still some works to 
do for the model (11. 3p . The first one is that, what is the asymptotic spreading 
speed when spreading happens? Since there is no other choice except spreading 
and vanishing, the second one is that we want to know the necessary condition for 
the disease to spread or to vanish. 
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